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After recalling a general non-perturbative expression for the luminosity-redshift relation holding in 
a recently proposed "geodesic light-cone" gauge, we show how it can be transformed to phenomeno- 
logically more convenient gauges in which cosmological perturbation theory is better understood. 
We present, in particular, the complete result on the luminosity-redshift relation in the Poisson gauge 
up to second order for a fairly generic perturbed cosmology, assuming that appreciable vector and 
tensor perturbations are only generated at second order. This relation provides a basic ingredient 
for the computation of the effects of stochastic inhomogeneities on precision dark-energy cosmology 
whose results we have anticipated in a recent letter. More generally, it can be used in connection 
with any physical information carried by light-like signals traveling along our past light-cone. 

PACS numbers: 98.80-k, 95.36.+X, 98.80.Es 

I. INTRODUCTION AND OUTLINE 

In a recent letter [1] (see also [2]) we have computed the effects of a stochastic background of inhomogeneities on 
the dct.criuinat.ion of dark-energy parameters in precision cosmology. The outcome of that analysis lias been that 
such perturbations cannot simulate a substantial fraction of dark energy: indeed, their contribution to the averaged 
flux-redshift relation is both too small (especially at large redshift) and has the wrong z-dcpendcncc. Nonetheless, 
stochastic fluctuations add a new and relatively important dispersion with respect to the prediction of the homo- 
geneous and isotropic Fricdmann-Lcmaitre-Robertson- Walker (FLRW) cosmology. This dispersion is independent of 
the experimental apparatus, the observational procedure, or the dispersion in absolute luminosity. Given the present 
(and probably near-future) limited statistics of Supernovae data, this, together with other phenomena, may prevent 
a determination of £l\(z) down to the percent level using the luminosity-redshift relation alone. 

In [1] we have presented the main ideas and most significant results of the calculation which, essentially, proceeds 
in two successive steps. The first one is the computation of the luminosity redshift relation dr / (z) (or, cquivalently, of 
the flux ~ d^ 2 ) at second order in perturbation theory. The method used in [1], being gauge invariant [3], allows us 
to express the result in a convenient gauge in which perturbations are known to second order, the so-called Poisson 
gauge (PG) [4] 1 . The second step consists of performing the relevant light-cone/ensemble averages, as in [2], and in 
inserting a realistic power spectrum of stochastic perturbations. This gives their effect on dark-energy parameters at 
the quantitative level. 

In this paper we will present full details about the first stage of this two-step process leaving the details about the 
second one to a future publication [9]. One reason for doing so is that the calculation of dh is independent of the rest 
of the calculation, has an interest of its own (i.e. irrespectively of its subsequent application to light-cone/ensemblc 
averaging) and could possibly find many other applications in precision cosmology. Furthermore, the result presented 
here for di, is valid in general, i.e. for any given background model 2 . 

The paper is organized as follows. In Section 2 we specify the PG up to second order in perturbation theory. We 
then recall the definition and special properties of an adapted system of coordinates introduced in [3] dubbed the 
geodesic light-cone (GLC) gauge. We also give the connection between the two gauges up to second order. In Section 



1 Following the pioneering work of [5], has been already computed to first order in the longitudinal gauge (for a CDM model in [6], 
CDM and ACDM in [7] and for a generic model in [2]), and to second order in the synchronous gauge, but only for a dust-dominated 
Universe, in [8], 

2 Except if caustics form. It has been argued [10] that the area distance is modified when caustics are present inside the past light-cone. 
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3 we present the actual calculation of <1l and express the final result, in compact form, in terms of perturbations 
in the PG, of the observed redshift, and of the observer's angular coordinates. The long, full expression of can 
be found in the Appendix together with a recollection of our definitions. In Section 4 we first offer some physical 
interpretation of the various terms appearing in and then show how the final result can be averaged over the 
observer's past light-cone reproducing the formulae used in [1]. In Section 5 we summarize the results and draw some 
short conclusions. 

We note that, after we submitted our short paper [1] - and while preparing this one - another group [11] has 
submitted a summary of their own calculation of dh in the PG and for a ACDM model. Since their calculation and 
ours are very different (and obviously completely independent) comparing the final outcomes, for such particular case 
of a ACDM model, will provide a very useful test of this highly non-trivial, long and somewhat tricky calculation. 



II. FROM THE POISSON TO THE GEODESIC LIGHT-CONE GAUGE AT SECOND ORDER 



A. The Poisson gauge 



Let us consider a non-homogeneous space-time approximated by a spatially flat FLRW Universe plus scalar, vector 
and tensor perturbations. In the so-called Poisson gauge (PG) ([4]), a generalization of the Newtonian (or longitudinal) 
gauge beyond first order, the corresponding metric takes the following standard form in cartesian coordinates: 

ds 2 PG = a 2 (r)) (-(1 + 2<$>)dif + 2uj l dr)dx l + [(1 - 2*)% + %] efecW) , (2.1) 

where $ and W are scalar perturbations, is a transverse vector (d l uii = 0) and hij is a transverse and traceless 
tensor (d l hij = = h\). This metric depends on six arbitrary functions, hence it is completely gauge fixed. Up to 
second order the (generalized) Bardeen potentials $ and VP are defined as follows: 

$ = V> + ^ (2) , * = V ,+ i^ 2 ) , (2.2) 

where we have assumed no anisotropic stress in order to set VP = $ = ip at first order. In this paper we shall consider 
oji and hij as second order quantities, the idea being that, in inflationary cosmology, first order scalar perturbations 
dominate over the others for small slow-roll parameters. On the other hand vector and tensor perturbations are 
automatically generated from scalar perturbations at second order (see e.g. [12, 13]). 



B. The geodesic light-cone gauge 



For problems associated with the observation of light sources lying on the past light-cone of a given observer, 
it is convenient to identify the null hypersurfaccs on which the photons reach the observer with those on which a 
null coordinate takes constant values. For this reason we have introduced in [3] an adapted system of coordinates - 
defining what we have called a "geodesic light-cone" (GLC) gauge - in which several quantities greatly simplify [3] 
while keeping all the required degrees of freedom for applications to general geometries. 

Let's recall [3] that the coordinates = (r, to, 6 a ) (with a = 1, 2, 6 1 = 6, 8 2 = <fi) specifying the metric in the GLC 
gauge correspond to a complete gauge fixing of the so-called observational coordinates, defined in [14-16]. The GLC 
metric too depends on six arbitrary functions (T, a two-dimensional "vector" U a and a symmetric matrix jab), and 
its line element takes the form 

ds 2 GLC = T 2 dw 2 -2TdwdT + j ab (de a ~U a dw){de b -U b dw) . (2.3) 

In matrix form, the metric and its inverse read: 

/ -1 -r- 1 -U b /T\ 

9g V lc = -T- 1 , (2.4) 

\-(U a ) T /T T 7 ab / 

where = (0,0), E/& = (Ui,U2), while the 2x2 matrices j a b and r y ab = (7 Q b) _1 lower and raise the two-dimensional 
indices. Clearly to is a null coordinate (i.e. d^wd^w = 0), and a past light-cone hypersurface is specified by the 
condition w = constant. We can also easily check that d^r defines a geodesic flow, i.e. that {d v r) V v (<9 m t) = (as a 
consequence of the relation g TT = — 1). 
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In the limiting case of a spatially flat homogeneous FLRW geometry, with scale factor a, cosmic time t, and 
conformal time parameter rj such that dr) = dt/a, the transformations to the GLC coordinates and the meaning of 
the new metric components are easily found as follows [3]: 

t = t , w = r + 77 , T = a(t) , 

U a = , Jabd9 a d9 b = a 2 (t)r 2 {dd 2 + sin 2 6d4> 2 ) = ^ b LRW de a d9 b . (2.5) 

Even though we will be mainly using the GLC gauge for a perturbed FLRW metric in the PG, it is important to 
stress that it is always possible to choose the GLC coordinates in such a way that r and t of the synchronous gauge 
are identified like in the above homogeneous FLRW limit [2]. As a consequence we can easily introduce with t a 
family of geodetic reference observers which exactly coincide with the static ones of the synchronous gauge. We also 
remark that, in GLC coordinates, the null geodesies connecting sources and observer are characterized by the simple 
tangent vector = g^ v d v w = g^ w = — (5^T _1 , meaning that photons travel at constant w and 9 a . This makes the 
calculation of the redshift and of the area distance particularly easy in this gauge. 

Let us denote by the subscripts "o" and "s" , respectively, a quantity evaluated at the observer and source space- 
time position, and consider a light ray emitted by a static geodetic source lying at the intersection between the past 
light-cone of a static geodetic observer (defined by the equation w = w ) and the spatial hypersurface r = r s with r s 
taken momentarily as a constant. The light ray will be received by such static geodetic observer at r = t > r s . The 
redshift z s associated with this light ray is then given by [3]: 

( _ (k»U n ) s _ (d^wd^T). _ T(w ,T ,9 a ) 

{ S> (fc"u„) (d^wd^o TK,r s >) ' [ J 

We will denote by T,(w ,z s ) the two-dimensional surface (topologically a sphere) which lies on our past light-cone 
(w = w ) and corresponds to a fixed redshift (z — z s ). In terms of the r coordinate this will correspond to imposing 
the equation r = T s (8 a ,w a , z s ) enforcing (2.6). Hereafter t s will denote this (in general angle-dependent) quantity. 

As said, also the area distance dA, related to the luminosity distance d^ of a source at redshift z s by the Etherington 
(or reciprocity) relation [17] 

d A = (1 + z s y 2 d L , (2.7) 
takes a particularly simple form in the GLC gauge [3]. We begin recalling the definition of d A [18] 

^ §: • < 2 - 8 > 

where dfl a is the infinitesimal solid angle at the observer, and dS is the cross-sectional area element perpendicular to 
the light ray at the source. Let us then show that, in the GLC gauge, we have [3]: 

d\ = ^~ ■ (2.9) 
smO 

Indeed, ~f a b is nothing but the induced metric on the surface T,(w ,z s ) provided this is parametrized in terms of 
the two "world-sheet" coordinates £° = 9 a and otherwise given by w = w , r = T s (8 a , w , z s ). Using the standard 
definition of an induced metric: 

ind ® X>X ® x ( \ in 1 n\ 

manifestly independent of the spacetime coordinates being used, we simply find: 

lTb d =la b , (2.11) 

since w = w (independently of 9 a ) and the only non-zero entry for the metric with a lower index t is g TW . We 
can also argue that the area element computed on this surface is orthogonal to the null geodesies and can therefore 
be identified with the dS of (2.8). Indeed, consider the projection of the photon momentum along the constant- z s 
hypersurface (which in this gauge, by (2.6), corresponds to constant T): 
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and the particular linear combination 

nW = a d^w + j3 d^T , (2.13) 
defining the normal to £ lying on the same constant- z s hypersurfacc and thus satisfying: 

n<*->0"T = => (d T T) a = T(0 |t T0"T) . (2.14) 

(z ) 

One can easily verify that is exactly parallel to fc^ii. Finally, using their constancy along the null geodesies, we 
can also identify 9 a with the angular coordinates at the observer's position where, within an infinitesimal region, we 
can take the metric to be flat. Therefore, as promised, 

£ = -£=*M. = ^.. (2.15) 
dn d 2 6sm9 sin6» 

The above expressions for the area distance dA singles out the flux $ ~ d^ 2 = (1 + z s ) _4 ii^ 2 as an important, and 
extremely simple, observable to average over the 2-sphcrc Yi{w ,z s ) embedded in the light-cone: 

(d- L 2 )( Wo ,Z s ) = (l + , s )-4£^ = (1 + Zs) -4i_^ = (1 + Zsr 4 4. 

J db J dS A(w ,z s ) 

A(w ,z s ) = { rf 2 CV7 ■ (2.16) 

J~S(w ,z) 

Here A is the proper area of £(u> ,z s ) computed with the induced metric j a b 3 . Exactly the same result follows 
from the averaging prescription of [3], which uses an alternative (and equivalent) definition of the surface H(w ,z s ) 
through two constraints obeyed by the coordinates. Eq. (2.16) holds non-pcrturbatively for any space-time and was 
the starting point of the computation of the average flux presented in [1] (see Eq. (6) therein). It can also be written 
in an elegant form in which the flux of an inhomogeneous Universe is compared to that of a FLRW one: 

= (d F A LRW f ^ . (2.17) 

It is now straightforward to formally express the result in a different gauge by simply changing the GLC gauge 
coordinates into those of the chosen new gauge. In our specific case wc would like to express 7 a t in terms of PG 
perturbations and, for our physical application to dark energy, as a function of the observer's angles d% — 6®. Taking 
in (2.10) the coordinates and the metric to be those of the PG we have: 



where y^ = y^(w 07 T s (w , z s ,9 a ),9 a ) define now the surface £(u> ,z s ) in the PG coordinates y 11 . Unfortunately, it 
is not trivial to find the explicit form of the above relation. We have found the easiest procedure to consist of: i) 
expressing the GLC gauge coordinates in terms of the PG coordinates and metric; ii) imposing the condition that 
the sources lie on E; iii) finally, inverting the second order transformation to express the outcome in terms of the 9 a 
angles. The first of this three-step process, which is of general interest, is carried out in the following two subsections. 



C. The second order transformation for scalar perturbations 



We now generalize to second order the transformation between the GLC gauge and the PG already obtained to first 
order in [2] . Before carrying on, let us mention that using the GLC approach means that we are taking many physical 
effects into account already at the level of the metric. In this approach, the geodesic equations for the observer 
and light rays are solved non-perturbatively, and their solutions are expressed in terms of the r, w coordinates. The 
outcome is that physical phenomena such as redshift perturbations (RP), rcdshift space distortions (RSD), Sachs- 
Wolfe effect (SW), integrated Sachs- Wolfe effect (ISW), peculiar velocities, lensing and others, are manifestly encoded 



3 As well known from the Nambu-Goto action in string theory! 
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in the metric, and are derived from a coordinate transformation. This is different from the usual approach, which hrst 
takes some perturbed metric, and then solves the geodesic equations order by order to construct physical observables 
(see, for example, [6, 20]). 

Since, by our assumption, vector and tensor modes appear only at second order and, as a consequence, will be 
decoupled from scalar perturbations, we can neglect them momentarily and add them at the end. 

Considering only scalar perturbations and using spherical coordinates (r,9 a ) = (r, #,</>), the PG metric defined in 
Eq.(2.1) can be rewritten as 

g^ G ^a{r 1 )- 2 diag{-l + 2^.l + 2^,{l + 2^) 1 ^) , (2.19) 

where 7q & = diag (r -2 , r~ 2 sin -2 9) , $ = ip + \4>^ — 2ip 2 and 4* = ip + ^ip^ + 2ip 2 . Following [2] we compute the 
GLC gauge (inverse) metric through: 

9cLc(*) = ^^ P c(y) ■ (2-20) 

where, as before, we indicate with = (r],r,8 a ) the PG coordinates and with x" = (t,w,9°") the GLC ones. Let us 
also introduce the useful (zeroth-order) light-cone variables rj± = rj ± r, with corresponding partial derivatives: 



d n = d + + d- 



d r = d + — 8- 



d 



1 



(2-21) 



Using these variables we solve the four differential equations obtained from Eq. (2.20) for the components g T Q LC = 
— 1, goi,c = 0> 9glc = 0' by imposing the boundary conditions that i) the transformation is non singular around 
r = 0, and ii) that the two-dimensional spatial sections r = const are locally parametrized at the observer's position 
by standard spherical coordinates. 

To this purpose we also introduce the following auxiliary quantities: 

P(r,,r,e a )= ^ d V '^li;(r 1 , ,r,e a ) , Q( V+ , r?_, 6 a ) = f ' dx x, 8 a ) , (2.22) 

where, hereafter, we use a hat to denote a quantity expressed in terms of (77+, # a ) variables, for instance 
^(r,+ ,r 1 .,9 a )=iP(r l ,r,9 a ). 

The sought for transformation can then be written, to second order in perturbation theory and with self-explanatory 
notations, as follows: 



,(0) 



rv 



(i) 



.(2) 



d V 'a(r 1 '))+a(r 1 )P(r ] ,r,e a ) 



dr) 



0(2) _ ^2 + {drP) 2 + 7 a6 Q a p Q b p\ (?? / ) p) 0a) ^ (2 23) 



ga(0) _|_ ga(l) _|_ Q~a(2) 



dx 



'/+ 



V 3(2) + 0(2) + ^Q +Q + ^ab Q aQ ^q] ^ ea) 



dx [% b d b Q] (v+,x,6 a ) 



dx 



7 ac Cc + v> e + a° 



(r,+,x,8 a ) 



(2.24) 



(2.25) 



where rji n represents an early enough time when the perturbation (or better the integrand) was negligible. In other 
words, all the relevant integrals (i.e. for all scales of interest) are insensitive to the actual value of rji n . Furthermore, 
we have used the following shorthand notations: 



Ccfa+,M°) 

e°(»?+.M o ) 



ld c w^( V+ ,x,9 a ) = l I dud c ^ +^ d+Q + f f d e Qd f Q (v+,u,9 a ) , (2.26) 



<M a(1) ( v+ , x,9 a ) + 2dJ a ^ (r] + ,x, 9 a ) 







du [r o c d c Q]{v+,u,9 a ) + [ 7o ac c> c Q](r7+,x,^) , 



(2.27) 



'/+ 



dJPM (v+,x, 9 a ) (d d 9 a W Xj e a ) - 6°id + Q(r} + ,x, 9 a ) 

\llo c d c Q]( V+ ,x,9 a ) ( jf du d d [T e d e Q] {r, + ,u,9 a )\ - \ [d + Q % b d b Q] ( V+ ,x,9 a ) . (2.28) 
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Let us now compute the various non-trivial entries of the GLC metric. Using T 1 = —d^w d v T g^ G , we obtain 

?- 1 = J-(l + eW+e<V) . (2.29) 

In terms of quantities implicitly defined in Eqs. (2.23)- (2.24) we find: 

eW = d + Q-d r P , (2.30) 

= d n w {2) + -d,.)r (2) - ipd n Q - <t> {2) +2ip 2 - d r Pd r Q - 2ipd r P -^ h d a Pd b Q . (2.31) 

The full explicit expression for e^ 2 ' can be then written as follows: 

e (2) = 1 ^(2) _ (2)^ + f_ + I( 5j .p)2 _ + Q +Q) . Q r p + hyab {2 Q a p . Q b p + g aQ . Q bQ _ idaQ . Q bP) 
^(2) + p) + ^ Q +Q + ^a b g aQ . dbQ \ Qa) 



n- 



ckr 



1 r d „. °w) dr 



4> (2) - V' 2 + (drPf + lZ b d a P ■ d b p] (rf, r, 9 a ) , 



(2.32) 



where the variables (r/,r,9 a ) have been omitted for the sake of conciseness. The computation of the GLC functions 
U a gives: 



u a = - \-dJ a ^ + - 7o ab <9 6 T (1) - dj a ^ + -^d^ + -d r r^d r e a ^ 

y a a a 



(2.33) 



where -r^ 1 ^ 2 ), fif ' 1 )^ 2 ) are implicitly defined in the coordinate transformations (2.23), and (2.25). U a is a measure 
of anisotropy of space-time in GLC coordinates. Substituting in (2.33) the explicit values of and fiK 1 )^ 2 ), we 

obtain the following explicit expressions for U a : 



U" 



7o ab ( - d b P ) + d+ [ - I dx [ 7o ab d b Q] (n+,x, e a ) 



+ I - (V + 8+Q) d+[- dx [% b d b Q] (r,+,x, 9 a ) + (-</> + 25 r P - d+Q) - 7o ab d b 



17+ 



1 



2^7o 6 d b P - - 7o cd d c P / 9 d [ 7o ah d b Q] Q ) + (9+Q - d r P)7 a& d b P 



1 



o7o b / drf-Y+dt (2) - ^ 2 + (d r P) 2 + % d d c Pd d P (V, r, r ) 



+ (d+ + d_) 



0(77) 



rf.Z' 



7o ac Cc + v r + a c 



(2.34) 



Finally, starting from 7 ah = flfp^y), we find: 



i(f7) 2 7°* = 7^ (1 + 2V») + [7S c d c 9 b{1) + (a o 6)1 + 7^ (V> (2) + 4^ 2 ) - d v 9 a ^d v 9 b{ ^ + d r 8 a{1) d r 9 b{1) 



2i/j 



fS c dJ b{1) + (a^b) + ^ d dJ a{1) d d 6 b{1) + lS c d c 9 b{2) + (a ^ b) 



(2.35) 
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More explicitly, in terms of the quantities defined in (2.22) and in (2.26)-(2.28): 



<il?l ah = 7o b (1 + 2V) + \ \lf jf ~ dx d d [f c d c Q] ( V+ ,x, 9 a ) + (a^b) 

\ rn- 
'v+ 



+ (i^ + 4^ 2 ) 7o b - \ lo c dcQ d + i- " dx [ 7o M d d Q] ( V+ ,x, 9 a ) + (a o b) 



+ V { 7o d / dx d d [j b c d c Q] {V+ ,x,6 a ) + (a^b) 

rv- 



dx d c [% e d e Q] 



7o~ / d c 

v+ 



( V+ ,x,9 a ) + (a^b) 



(2.36) 



D. The second order transformation for vector and tensor perturbations 



As already mentioned we can add the contributions of the tensor and vector perturbations by considering them 
separately. Using spherical coordinates (r,9 a ) — (r,0,4>), the tensor and vector part of the PG metric defined in 
Eq.(2.1) can be rewritten as: 



ds 2 PG = a 2 (rf) [-dr] 2 + 2v i dr)dx i + [("f ) tJ + x<j]da;W] , 



(2.37) 



corresponding to: 



(2.38) 



where rfo = diag(l,r~ 2 ,r~ 2 (sin 9)~ 2 ) is the (inverse) flat 3-metric. Here v 1 and x^ are the vector and tensor pertur- 
bation in spherical coordinates equivalent to the more standard definition and h lJ used in cartesian coordinates. 
They satisfy V^ 1 = ViX u = and (7o)yX^ = with V< the flat covariant derivative in spherical coordinates. 
Proceeding as for the scalar part of the metric, we first note that r is not affected by vector and tensor perturbations, 



while the light-cone coordinate w is: 



v 

dr/a(?/) , 

'»7»n 



w = r, + +Q^(r, + ,r,-,e a ) 



(2.39) 



(2.40) 



where 



Q {a) (r] + ,j^,9 a ) = dx a r (r] + ,x,9 a ) with a r 



v X 
2 4 



(2.41) 



As before, hats mean that (77+, rj-, # a )-coordinates are used. Finally, we find: 



r = frW+frW = 9 a + dx U a (r, + ,x,9 a ) - X ra (v + ,x,9 a ) + *ff(r, + ,x,6 a ) J* dy d b a r ( V+ , y, 6 a ) j . (2.42) 
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Following the same steps as for scalar perturbations, we then compute the non-trivial entries of the GLC metric: 



«(l)T-' = 1 - -v r - -x rr + d+ J dx <%+,*, a ) 
a(v) 2 l ab = 7o b - X ab 



(2.43) 
(2.44) 



^ jT~ dx d c U b (r 1+ ,x,e a ) -x rb (v + ,x,O a ) +*f b o d (v + ,x,6 a )J* dy d d a r { n+ ,y,e a )\ + (a o 6) 



C/ a = -<u a + a, ; a(2) = -v a + dJ a{ - 2) + <9 + (? a(2) = i ( -w a - X ra + lo b dx d b a r (r) + ,x,e a ) 



7+ 



rf.2' 



^(7 ?+ ,x,0 Q )-x ra ^+,a-^ a )+7o ,, (r/ + ,a ; ,^) / fta'^, y, a ) 



(2.45) 



where, again, the variables (rj, r, 9 a ) have been omitted for the sake of conciseness. Of course, these vector and tensor 
corrections to the FLRW metric have to be added to the scalar ones of the previous subsection. 



III. DETAILED EXPRESSION FOR d L (z,8 a ) 
A. The scalar contribution 



We now apply the above coordinate transformations to find the final expression of the luminosity distance in terms 
of perturbations in the PG, of the observed redshift, and of the observer's angular coordinates beginning, once more, 
with the scalar contribution. From Eqs.(2.7) and (2.9) we have: 



d L = (1 + z s ) 2 7 1 / 4 (shitf)" 1 / 2 . 
Let us start with 7. For a source emitting light at time rj s and radial distance r s we obtain from Eq. (2.35): 

7- 1 = det 7 afc = (a 2 r 2 shitf)" 2 {l + H's + 2<9 a (9 a(1) + 2^ + 12V' S 2 + 2d a 8 a{2) - A l0ab d+9 a{1) dJ b{l) 

+ %i) s d a 6 a{1) + 2<9 a a(1) <9 fc b(1) - dJWdbfPW } . 

We also need the expression for sin# up to second order in perturbation theory. This is easily given as: 



smv = suit 



Using Eqs.(3.2) and (3.3), Eq.(3.1) yields: 



(1 + z s ) A (a s r s ) { 1 -ip s - J 2 - \i> 



(2) 



- k 2 + ^ s j 2 + h j 2 ) 2 + — l — 

I 2 A sin 



where: 



J2 =2 



COt flW + dj" (1) J = -V o 0° (1) , K 2 = - 



cotee^ + d a 9 a ^ = ^v a <H 2 > 



(3.1) 



(3-2) 



(3.3) 



(3.4) 



(3.5) 



All the above quantities are evaluated at the source (apart from ip s , we neglect the suffix s for simplicity). 

At this point, for the explicit expression of the luminosity distance d^ at constant redshift, we need the first and 
second-order expansion of the factor a s r s = a(rj s )r s appearing in Eq. (3.4). To this purpose, we start from the explicit 
expression for the redshift parameter z s (see Eq.(2.6)), considered as a constant parameter localizing, together with 
the w = w condition, the source on H(w ,z s ). We then look for approximate solutions for rj s = rj s {z Sl 9 a ) and 
r s =r s (z s ,O a ). 
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Let us first define the zero-order solution rji°^ through the (exact) relation 

rf) 1 



l + z s 



(3.6) 



where a Q = a(?] ). Inserting now the result (2.29) into Eq. (2.6) and expanding a(rj s ) and T 1 with respect to the 



background solutions and r^' (where we define rj s = rf^' + r/g 1 ' + r/f and r s = 



,(0) 



_ J°) _i_ JX) , „(2) 



_ JO) , (1) . (2) 



ri + r 



r s ), we obtain: 



1 + Z S 



a(r) ) 



1K2 



where H s ~ a ^ ( s , - and 

o(»7. ) 



(3.7) 



^W^OX) , e W=eW(r ? (°),rf,^) , e?> = e^O,* - ) , 4 2) = e^^.r™/) , (3. 



,(l->2) _ 



Similarly, in order to compute r s and ri 2 \ we need to expand the w = w constraint by writing: 



where 



Wo = {r^ + r«> } + {r^ + r« + «,« } + + r™ + w& + } 



(3.9) 



(3.10) 



(3.11) 



,(l->2) 



(3.12) 



The additional terms ei 1 ^ 2 ' , tUj 1- '' 1 ' appearing in the above equations stand for the second order contributions coming 



from Taylor expanding ei X \ Ws , around the background source position 4 . More precisely, they originate from the 
fact that, at first order, quantities that are already first order are integrated along the unperturbed line of sight, while, 
at second order, first order terms have to be integrated along the perturbed line of sight. 
From Eq.(2.30) we obtain 5 : 



with 



^ - e« = J 



J = ([d+Q] a -[d+Q] )-([d r P] s -[d r P] o ) 



(3.13) 



(3.14) 



and where, for example, the term [d r P] s denotes the expression of d r P with r\ and r replaced by and ri°\ we 
also remark that [<9+Q] Q = —ij) . Then, from Eqs.(3.7), (3.10), recalling (2.24), we compute: 



(i) 



J 



-(0) 



J 



(3.15) 



These expressions are in accordance with our previous work [2]. We wish to note, already at this point, that the e 
terms correspond to redshift perturbations (RP). The first order term, e^, gives rise to the Doppler effect due to the 



4 There is no equivalent contribution at the observer position as r/o and r = are fixed quantities, with no perturbative corrections. 

5 In this paper we use Q s instead of the quantity ty a v introduced in [1], the two are directly related by Q B = — 2Ar]^ av . 
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peculiar velocities d r P, and the SW and ISW effects are combined together in d+Q s ~ d+Q a . Let us, in fact, recall 
that d r P can be rewritten as [2] 

d r P = v-h, (3.16) 
where n is the unit tangent vector along the null geodesic connecting source and observer, and where 



(3.17) 



are the "peculiar velocities" associated to a geodesic configuration perturbed up to first order in the PG. 

Let us now move to the second order quantities appearing in (3.8, 3.9) and (3.11, 3.12). We simply have, from 
Eq.(2.32), 



e (2) _ e (2) 



\ (4 2) - ^) + \ (4 2) - V4 2) ) + \ (tf - €) + \{[drP]s) 2 \{[drP]o? - + [d + Q] s ) ■ [d r P} 6 

4(7o 6 ) S (2d a Ps ■ d b P s + d a Q s ■ d b Q s - 4d a Q s ■ d b P s ) - \ lim [^ b d a P ■ d b P] 



1 

2 

1 

2 

1 

4 ./r,<°» + 



a(r]o) 



dx d-L 



4> (2) - i> 2 + (d r pf + js b d a p ■ d b p\ (?/, A^, e a ) 
4tn -i> 2 + (d r pf + ^ b d a p ■ d b p] {rf, o, e a ) 



p) + ^(2) + ^Q +Q + ffi> . . g bQ \ (j? (0) +) Xj 0») 



while 



We then have: 



i^ 2) = Qs {-[dlQ] s + [d+$] s + [d 2 P] s ] + -L {[d n i>] s + H s [d r P] s + [d 2 P] s } . (3.19) 



(3.18) 



(l-t-2) 



H 



-*jj s J + Q s (fa- [d+Q] s ) 



Using (3.7), (3.10), and recalling (2.24), we can now calculate rj s and r s to second order obtaining: 



vi 2) = 



~ (Hi + n 2 s )^ + & - - + ( £ W) 2 - e^e^ + H^H^ - e CD) 



H 



i (h'+h 2 j 2 



Us m 



and 



,.(2, 



+ (-v> s + [a+Q] a ) Q s 



(3.20) 



(3.21) 



,{0)- 



4 /„«>)+ 



dx 



£(2) + ^(2) + 4?/ 3 9+ g + -aft ^gl (^(0)+^ 



(3.22) 



In the second order terms we have the expected couplings between first order terms as well as the (also expected) 
genuine second order SW and ISW effects such as (ipi 2 ^ — ipo) and J dxd+tp^ . However, at second order, new effects 
come into play: most notably the tangential peculiar velocity d a P, the tangential variation of the photon path d a Q, 
and a RSD due to the peculiar acceleration d 2 P. The somewhat surprising appearance of tangential derivatives in rji 2 ^ 

(2) 

and r s is simply a reflection of working on a fixed- z surface. As a consequence, rcdshift perturbations originating 
from those of r, cq. (2.23), feed back on i] s ,r s and, eventually, on c?l(z) . 
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To conclude, combining these results, we obtain: 



a{ris)r s 
a(^ 0) )A77 



= 1 



At] 



1 - 



HI 



H s Ar, 



lp s J 



+ E s (^o + [drP]o) J + (-i>o - ^ + [d r P]s - \d r P]o) |^ 



4Ar? 



dx 



^(2) + ^(2) + 47/ 3 a+ g + ( ?? (0) + ) Xj 0a) 



(3.23) 



where 



(3.24) 



Let us also note that in Eq.(3.4) there are two other first order terms that have to be Taylor expanded up to second 
order around the background solution connected to the observed redshift z s , i.e. ip s and J 2 - We find: 

rPs=ti 1) +ti 1 ^ 2) =V'(^°\Ar ? ,n + ^[^-^](^ ) ) Ar ? ,n-Q S M(^ 0) J A J 7,^) , (3.25) 

Tin 



where we have used the 2-dimcnsional Laplacian A2 = dg + cot 9 dg + (sin0)~ 2 <9|. 

Collecting all the results obtained up to now, and inserting them in Eq.(3.4), we write our final result on the effect 
of scalar perturbations in the following concise form: 

d L (z s ,6°) d L (z.,0°) = 1 + S W iZs>e a ) + § (2) iZs>n (3 27) 



(l + z s )a Arj d F L LRW ( 
where: 



4 2) (z s , e a ) = [E S J -gj (V>« + J 2 (1) ) - 4^ 2) 4^ 2) + X {2) + . (3.28) 

Here ^\ ip^~* 2) , 4^ an d J^ 2) are implicitly defined in (3.25), (3.26) and and are the second order 
terms appearing in (3.4) and (3.23), namely: 

x (2) = -Ui 2) - Ul -k 2 + m + \{J 2 ) 2 + — l — a (e^) 2 + ( 70 ) ab d + 6 a Md-9 b u + laJWdtPM , 

2 2 2 4 sin 9 V ' 4 

Y<» = E s (e[V - 42) + e (i-2)) _ _2_^ sJ + Es + [drP]o) j + {[ g rP]s _ ^ o _ ^ _ [drP]o ) |l 



U s Ar] V 'H 2 / 2 4Ar? ^(0)+ 



0(2) + ^(2) d+Q + 7 afc 9 a Q (r ? (°)+,x,0 a ) . (3.29) 



Let us briefly point out that in dh several terms look similar to the ones that affect the shear at second order. In 
particular, following [19], the standard Born correction and lens-lens coupling are similar to the terms present in Jf, 
(7oU<M a(1 ><9_# 6(1 > and d a 9 b ^d b 8 a ^\ 

On the other hand, as already stressed, photons reach the observer traveling at constant 6 a . Therefore, the observer's 
angles are given by the 9 a which coincide with 6 a at the observer position but not at the source, hence dh should be 
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written in terms of 9 a rather than of 9 a . As a consequence let us consider the inverse form of Eq.(2.25): 

1 rv- 



ja(0) + ga(l) + ga{2) 



'/+ 



dx^\ V+ ,x,9 a ) / dy d b iP(r, + ,y,6 a ) 



dxr o b (v+^,0 a ) / dy d b iP(r, + ,y,6 a ) 



dx 7 ad (r/ + , x, ¥) / dy d d ^(r, + ,y, 6 a ) 



dx 



7o c Cc + «" + A' 



(3.30) 



The luminosity distance dL(z s ,9 a ) will then be given by Taylor expanding dL{z s ,9 a ) around 9 a (wc use a bar to 
denote that the luminosity distance is now expressed in terms of 9 a ). Using Eq.(3.30) we obtain: 

d L (z s ,9 a ) _ d L (z s ,9 a ) 2a ) + m ( §a) 

(l + z s )a Ar,- df L « w (z s ) ~ +6s ( s ' )+i>s { s ' ' ' 



with s ( s\z a ,e a ) = 4 1) (z s ,0 a ) 



S^(z s ,9 a ) = 6%\z s ,9 a )+d b \s^(z s ,e a ) 



(3.31) 



Equations (3.28, 3.31), supplemented with the vector and tensor contribution discussed in the next subsection, are 
our main result. More explicit expressions, where terms with different physical meaning are collected separately, are 
presented in the Appendix. 

B. The Vector and Tensor contribution 

Following the procedure just presented for scalar perturbations we start from the general expression for dh consid- 
ering now just vector and tensor perturbations 6 . We obtain: 



\l + \lhoUx ab Kv s ,r s ,9 a )- 4 Q) -i / 



4-/,+ 



where in terms of the quantity Q^ a > defined in (2.41): 



J. 



dx 



nt f 



(0) + 
T)i - X 



7 A 2 Q^(r,^ + ,x,9 a ) 



(3.33) 



and where a s r s is a quantity that still needs to be expanded with respect to the observed redshift. In order to do 
that, we first write the analog of (3.7): 



1 



1 + 2., 



1 + 



V Y 
h — 

2 4 



(v s ,r s ,e a ) 



v X 
2 4 



(r] o ,0,e a )-jM 



(3.34) 



with 



J 



(a) 



dxd + a r (r,t,x,e a ) = a r + d+Q^ 



(3.35) 



Expanding rj s as r/ s = + rjs and imposing (1 + z s )a(rjs V> ) = a 0l we get 

1 



v y 

h — 

2 4 



(vo,o,e a ) 



v X 
2 4 



{4 0) ,rW ,9 a ) + jW 



(3.36) 



Note that an expression for the contribution of vectors and tensors to has been derived recently in [20]. 
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Using also the transformation of w, Eq. (2.40), we get the expression of r s : 



At? ■ 



1 



v X 
2 4 



v X 
2 4 



fe,o,n-^ (a) \ -Q (a) 



(3.37) 



and finally reach the conclusion that the luminosity distance at linear order in vector and tensor perturbations 
(regarded themselves as second order quantities) is 



d (v,T) 

d FLRW 



i+4!t = i- 



Q 



Ar? 

(0) + 



T + l) (I + A) + J(a) ) + ^°Ux ab M°\r?\8«) 



{ lv a [v a -x ra ](vi 0)+ ,x,e a )+ f 



(^° )+ - x) 



rA 2 Q^(r ? (°)+,x,^) 



Using the transvcrsality and trace-free conditions on the perturbations: 



► (3.38) 



V a V a = - [d r + -)v r , VaX ra = -[dr + -)X 



(ioUx ab = -x r 



we finally get an expression that depends only on v r and x r 



d {V,T) 
d FLRW 



Q 



1 

4 /„(°)+ 



%,Xo_ 

2 4 



dx 



4 X S 



d r + - ) X T 



1 



A 2 g(°) 



, x, I 



(3.39) 



(3.40) 



where one should interpret r 



— — inside the last integral. 
We note, once more, the nature of the terms appearing in (3.40); in the first line we see a SW term as well as 
an average/integrated SW effect for the vector/tensor perturbation. The second line involves frame-dragging and a 
"magnification" term for tensors/vectors proportional to the laplacian of the perturbation on the 2-sphere. 
Our final expression for is thus: 



d L {z s ,9 a ) 
d F L LRW {z s ) 



1 + ~5^\z s ,6 a ) + 6%\z s ,9 a ) + 5 ( *> T (z s , 6 a ) 



(3.41) 



where we replaced 6 a with 6 a in Sy ^ to get Sy' T since this is considered already as a second-order quantity. 



5(2) 



IV. INTERPRETATION OF d L (z,6 a ) AND APPLICATION TO THE AVERAGED FLUX 

In the previous Section we have obtained a "local" expression for di(z,0 a ), expression that can find a number of 
possible applications. Note the importance of giving the result in a gauge which is convenient in terms of computing 
(or just writing) cosmological perturbations (here the PG) but also of expressing the final outcome in terms of the 
GLC angular coordinates, since, given the constancy of the 6 a along the null geodesies, these correspond to the 
observer's angular coordinates. 

In this section we will first make some comments on the physical meaning of the various terms appearing in our 
final result. Finally, we will make contact between the local expression of d^ and its angular and ensemble averages 
stressing how those of <i£ 2 (hence essentially of the flux) lead to the expressions used in [1]. 

The different terms appearing in 6g , \ Oyr can ^ G rou S nr y classified as follows: 

(2) (2) 

• Redshift Perturbations. These are the e s — e terms as well as the analogous vector/tensor v r , x rr terms 
appearing in (3.34). As mentioned in the text, at first order they include Doppler effect of peculiar velocities 
(d r P s — d r P ), SW and ISW (in agreement with the results obtained in [6]). At second order additional effects 
such as the tangential peculiar velocity d a P, the tangential variation of the photon path d a Q, and RSD also 
appear. Peculiar velocities are the dominant contribution at low redshift z < 0.2, when we average generic 
functions of the luminosity distance di,{z) (see [1, 2]). 
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• Perturbed trajectories. At second order, first order integrated quantities are evaluated along the perturbed 
geodesies giving rise to (1 — > 2) terms. 

• SW and ISW effects coming from the evaluation of the area distance. Once again, in our notation, they are 
simply combined as (d+Q s — d+Q ). There is also an equivalent effect in the tensor/vector contribution. 

• Lensing. These are the magnification J2, K2 terms as well as the shear in d a 9 b ^ db9 a ^ 1 \ They are the most 
important contributions at high redshifts z > 0.5, when we average generic functions of the luminosity distance 
d L (z) (see [1, 2]). 

• Frame dragging, in the vector contribution. 

Let us now consider a possible application of these results, application already presented in [1]. We first note that 
the vector and tensor contributions vanish when we average a function of over the angles. Indeed, if for each Fourier 
mode we choose our z-axis in the direction of the wave-vector, and we expand the vector and tensor contributions in 
spherical harmonics, their contribution turns out to be proportional to e ±l< ^ and e ±2l< ^, respectively. In both cases 
their angular integration will give zero. 

We have already seen in Section IIB how the averaged flux takes the very simple form of a fraction (Eq. (2.16)) 
where the numerator is simply a pure number (basically the observer's solid angle and, in the denominator, we 
have the invariant area of the £(iu , z s ) surface. In order to evaluate the latter in terms of the PG metric perturbations 
we will start expressing ^F) in the Poisson gauge while still using the angular GLC coordinates (as done in the previous 
section for dc). 

Starting from Eq.(3.2) we can obtain, in a straightforward way, the following expression: 

V7 = (a s r s ) 2 (sin0) jl - (2^ s + dj a ^) + [ - V4 2) - dj a ^ + 2( lQ ) ab d + 9 a ^d-9 b ^ + 2^ s dJ a ^ 

+ \d a 9 a{1) d b 9 b ^ + l0 o flWD&flMD] I , (4.1) 

where, in particular, a s r s is given by Eq.(3.23). Next we express sin# in terms of 9 a (angles seen by the observer). 
Starting from Eq.(3.30) we obtain: 



sin 9 = sin 9 



l + cot0(0«+0< a >)-i(0M) a 



(4.2) 



Similarly we can Taylor-expand the rest of the terms present in Eq.(4.1) around 9 a (using Eq.(3.30)) and around the 
background values rji " 1 and ri ' , and arrive at an explicit form for y^y. We omit writing the explicit - and not so 
illuminating - expression. This can be finally integrated over the 9 a angles, according to Eq. (2.16) with £ a = 9 a , to 
obtain A(w , z s ). The final result can then be put in the form of Eq. (8) of [1], namely: 

7,(, s ) = ( a (^))A ?? )- 2 ^^ = |^sin0(l + I 1+ I 1 , 1 +I 2 ) . (4.3) 

where one can easily show that the following connection should exist between the various quantities appearing on the 
r.h.s. of (4.3) and those in (3.31): 



II = 2S ( g> + (t. d.) (l) 

f ] + (4 i; 



l hl +l 2 = 2~5f + (4 X) ) 2 + (*• d -) (2) , (4-4) 



where the (t. d.)*- 1,2 ^ appearing in (4.4) denote total derivatives terms w.r.t. the 9 a angles giving vanishing contribution 
either by periodicity in <f> or by the vanishing of the integrand at 9 = 0, 7r. As an example of such terms consider the 
first order contribution Ti whose explicit expression is: 

I 1 = -2^V^n + 2(W-i-Q s ) . (4.5) 

This expression can be compared with the one of given in (3.28). Apart from an obvious factor two, the 

expression for X\ lacks the term which is precisely a typical one that vanishes upon angular integration. 
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Still dropping irrelevant total derivatives, the two second order terms appearing in (4.3) take the following explicit 
form: 



= 2S S I \ [tf - + \{[d r P] s ) 2 - \{{d r P] ) 2 - + [d+Q].) ■ [d r P] 6 
1, „ „ „ „ „ „ 1 



-( 7o ab ) s (2d a P s ■ d b P s + d a Q s ■ d b Q s - Ad a Q s ■ d b P s ) - - lim [^ b d a P ■ d b P] 

4 Z i — >0 

-[d 2 + Q] s + [d + ii] s + [d 2 P] : 



+ 



+ + H a [drP]. + [d 2 r P]> 

ris 



+ 



"2 



drf-^Qr-dr [~^ 2 + (d r P) 2 + 7 ab 9 Q P • d b P] (?/, At?, 6 a ) 
Vin a(rf s ') 

dn'^Qrdr [-V- 2 + (drPf + lfd a P ■ b P] (r/, 0, 0°) 

(0)- 



dx d+ 



HsAr) 



1 



V; d+Q + - 7o a6 • d a Q ■ d b Q 



K 
n 2 



J - ±i> s J + 2S S U° ~ + [d r P]o )J + ^ 



A?/ 

0, , 



+ 2 (V> s - + [d r P]s - [d r P]o) ^ + (70 )^aQ^ 6 ^ + 



n s 



1 1 g 

8 sin 90 



Ary 
) M 



2 K 



2 



rZ.i' 



cos 



>?i 0) - 

(0) + 



A?y 7 J) (°)+ 
dx {f b d b Q}( V ^+,x,d«) 



Q 
Ar) 



4> d + Q + - 7o a6 d a Q 8 b Q 



(4 0)+ ,xJ a ) 



(4.6) 




where it is important to stress that all these quantities have their angular dependence expressed in terms of 8 a . Let 
us also point out that the last term in Eq.(4.6) corresponds to a total derivative and thus to a boundary contribution 
that superficially looks non vanishing. We believe that this is the result of a naive treatment of the angular coordinate 
transformation which becomes singular near the poles of the 2-sphere. This contribution has indeed the same form 
as that of an overall 5*0(3) rotation connecting 9 a and 9 a . Modulo this subtlety, one can explicitly check (through a 
long but straightforward calculation) that Eq. (4.4) is indeed satisfied. 

To conclude, using the results (4.5-4.7) in (4.3) and considering the ensemble average (see, for example, [21-23]) 
of {dJ 2 ){w T z s ) for a stochastic spectrum of inhomogeneities, we obtain the results already discussed in [l] 7 . As 
anticipated, this last, more phenomenological step, will be described in detail in a future publication [9]. 



The observational consequence of the use of the ensemble average and of a stochastic spectrum of inhomogeneities were also recently 
considered, in a different context, in [24]. 
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V. CONCLUSIONS 

We have presented an explicit calculation of the luminosity distance as a function of the redshift and angular 
coordinates measured by a geodetic observer. The result, being expressed in terms of the Poisson-gauge metric 
perturbations up to second order, is a suitable starting point for determining the quantitative effects of cosmological 
perturbations once a particular inhomogeneous model is chosen. 

Our approach is making heavy use of a newly introduced [3] geodetic light-cone (GLC) gauge endowed with some 
characteristic and extremely useful properties. Indeed, both the redshift and the luminosity distance are simply 
expressible in terms of the GLC metric while the past light-cone of the observer reduces to fixing one (null) GLC 
coordinate. Furthermore, since the null geodesies going from the source to the observer are at constant GLC gauge 
angles 9 a these can be identified with the observer's angular coordinates with respect to which various moments can 
be in principle computed along the lines already discussed in [6, 20]. 

The advantages of the GLC gauge have been illustrated here for the case of averaging the flux $ ~ d^ 2 whose 
interest for the determination of dark-energy parameters has been already discussed in [1]. In this case the problem 
is essentially reduced to the evaluation of the proper area of the fixed z s surface lying on our past light cone. This 
simple result can be applied to fully deterministic (classical) inhomogeneous models (such as LTB models of the kind 
discussed in [25]) even when relaxing a fine-tuned condition on the position of the observer, or, more realistically, to 
the stochastic inhomogeneous models that follow from inflation as done in [1, 2]. 

Since the luminosity distance is related to the magnification of an image, our results could potentially also have 
consequences in studies of weak lensing surveys or on the "anti-lensing" effect due to a stochastic distribution of large 
voids [2 (>]. Another application could be to the determination of dark-energy parameters via the so-called redshift 
drift (see, for example, [27]) as already anticipated in [3], or to the analysis of CMB anisotropics, including non- 
gaussianity, B polarization due to tensor modes, etc. More generally, our approach can be useful whenever dealing 
with information carried by light-like signals travelling along our past light cone. 
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Appendix. Detailed expression of 8^ (z s , 9 a ) 

The second order corrections appearing in (3.31) can be conveniently grouped as follows: 

5f (z„ 9 a ) = 6% th + 8f} s + 4t ed (A.l) 
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where S^J th is for the terms concerning the photon path, 6p OS for the terms generated by the source and observer 

— /n\ 

positions, and SjJtxed * s a m i xni § or * both effects. Their explicit expressions are: 



<& = s.{ - \ - 0i 2) ) + \ (^ 2) - ^) + \ (v, - w 2 - <m 



(i) 



+ (ip -i> s - A 1] )[d+Q] s + ^(j$ b ) s d a Q s ■ d b Q s + Q s (-[d 2 + Q} s + [d+j,].) + ^-(ij + [d+Q] s )[d v tb] 6 



2 / d7? T^h r 



ii,„ 



i 



4 



da; 9+ U( 2 ' + ^ + 4</> 9+Q + % b d a Q d b Q\ 9 a ) 



,(0)- 



^ (l - |f ) ^(Vo + [d + Q] s ) 2 ^AzM^o + [d + Q]s) 

-d (lb +J {1) + 0± 



(0) + 



1 



>7< 0) - 

(0) + 



[7 ah d b Q] (vi 0)+ ,x,6 a ) 



16 \ .lrji 0) + 



dx [j bc d c Q](r 1 ^ + ,x,6 a ))d b [ /'' q S )+ dx [jf d d Q] ^J a ) 



1 



4A?y J^to) 



dx 



0(2) + V }(2) + 4 ^ a+ g + ^6 g a g ^gj (^(0)+ ja!j 0«) 
1 



+ wr(1>o + [d + Q] a )\-[W] a + [arii>].+ 



H 



+ Q s < [d r ip] s + d+ 



Arj 2 J v (o) 



dv'A 2 rl;(r ] , ,r ]o -r 1 , ,e a ) 



dx- 



(^ 0)+ - x) 2 J v i 0)+ 



dyA 2 ^r,^+,y,e a ) 



i r n 



2Ar] 2 J (o) 



dr)'A 2 i>(v',Vo-v',O a ) 



1 



16sin z 6> \-V c)+ 



[7 lfc d b Q](r,^ + ,x,6 a )) , 



(A.2) 



4- = (M- " I^o) 2 + hS b ) s d a Ps ■ d b P s - lim [jtd a P ■ d b P] - — {[d r P] s - [d r P] ) {U s [d r P] s + [d 2 r P] 

2 r^O tis 



dr] 



Hi 



d r [{d r Pf + ^ b 8 a P ■ d b P] (if, A Vl a ) + I dv'^Urd r [(d r P) 2 + ^ b d a P ■ d b p] ( v ',o, e a ) 



a(rj ) 
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f(2) 



2tp Q -ip s + d+i 



[d r P]o - ([drP]s - [drP]o) (^"[^]. " 4^ ) 



+ ^r(i>o + [d+Q}sMP} s + Q s ■ [d 2 r P] s + d a ([d r P} a - [d r P] ) ■ \ 



dx [% b d b o\ (T 1 i o)+ ,x,n 



1 

Ar] 



([d r P] s - [d r P] ) 



H s 



1 - 



K 

HI 



(ipo + [d+Q] s ) + —ip s + Q a 
rts 



+ y-([drP]s - [d r P]o) ■ { [d^]s - [drift, 



Arf J (o) 



d V 'A 2 i ) (r 1 ', Vo -r 1 , ,e a ) 



(A.4) 



The various quantities appearing in the above equations are defined in the main text but are reported again below 
for the reader's convenience: 



P(r,,r,6 a ) 



n- 



Q(r) + ,ri-,e a ) = / dx iP(r, + ,x,6 a ), 

Jrj + 

„ 1 



HsAr] 



J > ^ 2 



cot e §w + d a a{1 ^ = -V a 6i Q(1) 



cote § (2) + d a e a( - 2) 



4 1] = ^ / ... dr, 



i mo (0) 



A 2 ^( V , Vo - V ,9 a ), 



J = ([d+Ql - [d+QD - ([d r P] s - [d r P} ) with [d+Q} = -V^o 



I. Ben-Dayan, M. Gasperini, G. Marozzi, F. Nugier and G. Veneziano, [astro-ph/1207.1286]. 
I. Ben-Dayan, M. Gasperini, G. Marozzi, F. Nugier and G. Veneziano, JCAP 04, 036 (2012). 
M. Gasperini, G. Marozzi, F. Nugier and G. Veneziano, JCAP 07, 008 (2011). 

E. Bertschinger, Cosmological Dynamics in Cosmology and Large Scale Structure, proc. Les Houches Summer School, 
Session LX, ed. R. Schaeffer, J. Silk, M. Spiro and J. Zinn- Justin (Amsterdam:Elsevier Science), 273-347, 1996; 
L. Bombelli, W. E. Couch and R. J. Torrence, Class. Quant. Grav. 11, 139 (1994). 

M. Sasaki, Mon. Not. Roy. Astron. Soc. 228 (1987) 653; M. Kasai and M. Sasaki, Mod. Phys. Lett. A2 (1987) 727. 

C. Bonvin, R. Durrer and M. A. Gasparini, Phys. Rev. D 73, 023523 (2006) [Erratum- ibid. D 85, 029901 (2012)]. 

T. Pyne and M. Birkinshaw, Mon. Not. Roy. Astron. Soc. 348, 581 (2004). 

E. Barausse, S. Matarrese and A. Riotto, Phys. Rev. D 71, 063537 (2005). 

I. Ben-Dayan, M. Gasperini, G. Marozzi, F. Nugier and G. Veneziano, in preparation. 

G. F. R. Ellis, B. A. C. C. Bassett and P. K. S. Dunsby, Class. Quantum Grav. 15 (1998) 2345. 

0. Umeh, C. Clarkson and R. Maartens, arXiv: 1207.2109 [astro-ph.CO]. 

S. Matarrese, S. Mollerach and M. Bruni, Phys. Rev. D 58, 043504 (1998). 

N. Bartolo, S. Matarrese and A. Riotto, JCAP 0605, 010 (2006). 

R. Maartens, PhD Thesis, University of Cape Town, South Africa (1980). 

G. F. R. Ellis, S. D. Nel, R. Maartens, W. R. Stoeger and A. P. Whitman, Phys. Rep. 124, 315 (1985). 
C. Clarkson and R. Maartens, Class. Quantum Grav. 27, 12 (2010). 

1. M. H. Etherington, Phil. Mag. 15, 761 (1933). 

J. Kristian and R. K. Sachs, Ap. J. 143 (1966) 379. 



19 



[19] S. Dodelson, E. W. Kolb, S. Matarrese, A. Riotto and P. Zhang, Phys. Rev. D 72, 103004 (2005). 

[20] E. Di Dio and R. Durrer, Phys. Rev. D 86, 023510 (2012). 

[21] N. Li and D. J. Schwarz, Phys. Rev. D 78, 083531 (2008). 

[22] C. Clarkson, K. Ananda and J. Larena, Phys. Rev. D 80, 083525 (2009). 

[23] C. Clarkson and O. Umeh, Class. Quantum Grav. P28, 164010 (2011). 

[24] G. Marozzi and J-P. Uzan, arXiv: 1206.4887 (2012). 

[25] M. Celerier, Astron. Astrophys. 353, 62 (2000). 

[26] K. Bolejko, C. Clarkson, R. Maartens, D. Bacon, N. Meures and E. Beynon, arXiv: 1209.3142 [astro-ph.CO]. 
[27] C. Quercellini, L. Amendola, A. Balbi, P. Cabella and M. Quartin, arXiv:1011.2646 [astro-ph.CO]. 



